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Abstract. We consider the Polyakov loop-extended two flavor chiral quark-meson 
model and discuss critical phenomena related with the spontaneous breaking of the 
chiral symmetry. The model is explored beyond the mean-field approximation in the 
framework of the functional renormalisation group. We discuss properties of the net- 
quark number density fluctuations as well as their higher cumulants. We show that 
with the increasing net-quark number density, the higher order cumulants exhibit a 
strong sensitivity to the chiral crossover transition. We discuss their role as probes of 
the chiral phase transition in heavy-ion collisions at RHIC and LHC. 



1. Introduction 

Lattice Quantum Chromodynamics (LQCD) [1] has confirmed that at finite 
temperature and small baryon density QCD exhibits restoration of chiral symmetry 
and deconfinement. However, owing to the fermion sign problem, the thermodynamics 
of strongly interacting matter at large baryon densities is presently not accessible by 
first principle LQCD calculations. A viable framework for exploratory studies of QCD 
at finite density is however offered by phenomenological models and effective theories. 

The Polyakov loop extended Nambu-Jona-Lasinio (PNJL) and quark-meson 
(PQM) models reproduce essential features of QCD thermodynamics already in the 
mean-field approximation. However, to correctly account for the critical behavior and 
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scaling properties near the chiral phase transition, thermal and quantum fluctuations 
must be included in a non-perturbative manner. This can be achieved e.g. by using 
methods based on the functional renormalization group (FRG) [2]. 

We compute the cumulants of the net-quark number density iXn ) ^^'^ electric 
charge density [Xn) finite temperature and baryon chemical potential. We show 
that the higher-order cumulants exhibit a distinctive structure in the vicinity of the 
transition, where they become negative. Such a structure follows from the scaling 
functions of the three dimensional 0(4) universality class [31 IH |5]. Ratios of such 
cumulants have been suggested as sensitive probes for the critical behaviour and the 
rapid change of degrees of freedom at the QCD phase transition [6]. 



2. Polyakov loop-extended quark— meson model 



The quark-meson model is an effective realization of the low-energy sector of QCD, 
which incorporates chiral symmetry. By introducing the coupling of quarks to a 
uniform temporal color gauge field represented by the Polyakov loop, the model can be 
employed to explore properties related to color confinement. 

In order to account for non-perturbative mesonic fiuctuations in the PQM model, 
we employ methods based on the functional renormalization group (FRG). The FRG 
involves an infrared regularization of the fiuctuations at a sliding momentum scale 
/c, resulting in a scale dependent effective action F^, the so-called effective average 
action |2]. We treat the Polyakov loop as a background field, which is introduced self- 
consistently on the mean-field level while the fiuctuations of the quark and meson fields 
are accounted for by solving the FRG fiow equations. 

The fiow equation for the scale-dependent grand canonical potential density, VLk = 
TTk/V, for the quark and meson subsystem reads [71 [8] 



9,nfc(£,r;T,/i) 
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Here usiET^y^T) is the bosonic distribution function 

nB{E^y,T) 



exp(E^,^/T) 
E^ = 



where the primes denote derivatives with respect to p = (cr^ + 7f^)/2 oi Q = Q + ca. The 
fermion distribution functions N{i, t; T, p) and N{i, t; T, fi) = N{t, i; T, -p), 

1 -\- 2te'^^^''~^^ + ie^i^^^"'''^ 



N{iJ*;T,fi) 



(2) 



are modified because of the coupling to the gluon field. Finally, the quark energy is 
given by Eg = y/k'^ + 2g'^p. 
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The minimum of the thermodynamic potential is determined by the stationarity 
condition 



dQi 



da 



dfli 



da 



- c = 0. (3) 

The flow equation ([T]) is solved numerically with the initial cutoff A = 1.2 GeV (see 
details in Refs. O [8]). The initial conditions for the flow are chosen to reproduce 
vacuum properties: the physical pion mass = 138 MeV, the pion decay constant 
= 93 MeV, the sigma mass ra„ = 600 MeV and the constituent quark mass = 300 
MeV at the scale — )• 0. The symmetry breaking term, c = m'^f^^, corresponds to an 
external field and does not flow. 

By solving Eq. ([1]), one obtains the thermodynamic potential for the quark and 
meson subsystem, Qk^o^i, i*;T, as a function of the Polyakov loop variables i and 
i*. The full thermodynamic potential i*; T, /i) in the PQM model, including quark, 
meson, and gluon degrees of freedom is obtained by adding the effective gluon potential 
U{i,t) to l],^o(^,r;T,/i): 

Q{e,t;T,fi) = Qk^o{i,t;T,fi)+Uii,n. (4) 

At a given temperature and chemical potential, the Polyakov loop variables, i and i*, 
are then determined by the stationarity conditions: 

^^](£,^;^,/i) = o, ±.n{e,t;T,f,) = o. (5) 



3. Cumulants of conserved charge 

The fluctuations of conserved charges are quantified by cumulants Xn'*^; which are 
generalized susceptibilities obtained by taking derivatives of the pressure p = T'^p 
with respect to the corresponding chemical potential o = TfiBo: Xn''^ = fr^^ ■ 

The cumulants of conserved charges are sensitive probes of the chiral phase 
transition. They indicate the position, the order, and in case of the second-order phase 
transition the universality class of the corresponding phase transition. The net baryon 
number density is discontinuous at a first-order transition, whereas the susceptibihty 
X2 and higher cumulants diverge at the critical end point [9] and at the spinodal lines 
of a first-order chiral phase transition [lOj. In the chiral limit and at non-zero chemical 
potential, all generalized susceptibilities Xn with n > 2 diverge at the 0(4) chiral critical 
line inilTT], while at vanishing chemical potential this divergence shows up only for n > 4. 

Close to the phase transition the sixth and eighth order cumulants decrease from 
their positive values in the hadron resonance gas to negative values, in accordance with 
the expected dominance of the 0(4) singularity [3]. Since the fourth order cumulant is 
unaffected by chiral critical phenomena at zero chemical potential, it remains positive 
for small chemical potentials, as in the hadron resonance gas (HRG). However, at larger 
densities, it picks up a non-trivial contribution from higher order /i = cumulants, 
as seen in the Taylor expansion xf (/i) = xf (0) + |xf (0)/i^ + 0{fi'^). Thus, the fourth 
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Figure 1. Left: The chiral crossover line (dashed) and the first zeros in Xn for 
n = 4, 6 and 8. Middle: The parameter range for which the sixth order cumulant of 
net-baryon number fluctuations is negative. Right: The parameter range for which the 
sixth order cumulant of electric charge fluctuations is negative. 



order cumulant and consequently the kurtosis (xf/xf ) i^^y become negative close to 
the crossover transition even if the critical end point does not exist [3l [TTl [12] . The 
location of the first zeros of Xn n = 4, 6 and 8 are shown in Fig. 1. The middle 
panel shows that the region, where the sixth order cumulant of the net baryon number 
fluctuations is negative, is closely correlated with the crossover transition. 

The FRG method can also be applied to compute the fluctuations of electric 
charge in the PQM model. Here the differences relative to the mean field calculations 
are essential, since the FRG approach accounts for charged pion contributions to Xn- 
In Fig. 1 we show the range of negative xS ^ which is similar to that found for x% ■ Thus, 
negative fluctuations of the sixth order moments of net baryon as well as net electric 
charge fluctuations can both be attributed to chiral critical dynamics. Consequently, 
the experimental observation of negative sixth order cumulants in heavy ion collisions at 
RHIC and LHC would indicate that the chemical freeze-out takes place in the vicinity 
of the chiral crossover transition [3]. 

References 

[1] M. Cheng et al., Phys. Rev. D 77, 014511 (2008). 

[2] J. Berges, N. Tetradis and C. Wetterich, Phys. Rept. 363, 223 (2002). 

[3] B. Friman, F. Karsch, K. Redlich and V. Skokov, arXiv:1103.3511 [hep-ph]. 

[4] F. Karsch and K. Redlich, Phys. Lett. B 695,136 (2011). 

[5] F. Karsch, these proceedings |arXiv:1108. 08341 [hep-lat]]. 

[6] S. Ejiri, F. Karsch and K. Redlich, Phys. Lett. B 633, 275 (2006). 

[7] V. Skokov, B. Stokic, B. Friman and K. Redlich, Phys. Rev. C 82, 015206 (2010). 

[8] V. Skokov, B. Friman and K. Redlich, Phys. Rev. C 83, 054904 (2011). 

[9] M. A. Stephanov, Phys. Rev. Lett. 102, 032301 (2009) 

[10] C. Sasaki, B. Friman and K. Redlich, Phys. Rev. D 77, 034024 (2008). 

Phys. Rev. Lett. 99, 232301 (2007). Phys. Rev. D 75, 074013 (2007). 

[11] M. A. Stephanov, FarXiv: 1104. 16271 [hep-ph]. 

[12] M. Stephanov, these proceedings. 



